Density-functional study of Cu atoms, monolayers, and coadsorbates 

on polar ZnO surfaces 



CO 

o 
o 

CN 



o 

CZ3 



B. Meyer and D. Marx 
Lehrstuhl fur Theoretische Chemie, Ruhr-Universitat Bochum, 44780 Bochum, Germany 

(Dated: February 2, 2008) 

The structure and electronic properties of single Cu atoms, copper monolayers and thin copper 
films on the polar oxygen and zinc terminated surfaces of ZnO are studied using periodic density- 
functional calculations. We find that the binding energy of Cu atoms sensitively depends on how 
charge neutrality of the polar surfaces is achieved. Bonding is very strong if the surfaces are stabilized 
by an electronic mechanism which leads to partially filled surface bands. As soon as the surface 
bands are filled (either by partial Cu coverage, by coadsorbates, or by the formation of defects), 
the binding energy decreases significantly. In this case, values very similar to those found for 
nonpolar surfaces and for copper on finite ZnO clusters are obtained. Possible implications of these 
observations concerning the growth mode of copper on polar ZnO surfaces and their importance in 
catalysis are discussed. 

PACS numbers: 68.43.Bc, 68.43.Fg, 68.47.Gh, 82.65. +r 
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I. INTRODUCTION 

The ternary composite system Cu/ZnO/A^Oa is an 
important commercial catalysts for the methanol synthe- 
sis, the water-gas shift reaction and the methanol-steam 
reforming processi^ Although pure ZnO is a good high 
temperature catalyst for hydrogenation and dehydro- 
genation reactions, only the combination with Cu makes 
it an effective and selective catalyst for low temperature 
applications on an industrial scale. In particular the role 
of copper, including its electronic and atomistic structure 
on ZnO surfaces, is not clearly understood. There ex- 
ist many controversial models to explain how copper en- 
hances the activity of such ZnO-based catalystsi^ Dif- 
ferent "active species" are proposed in the literature, like 
the existence of Cu + impurities at the ZnO surfaces, but 
the most commonly accepted picture is that Cu forms 
small metallic clusters which are dispersed on defect-rich 
and heterogeneous ZnO surfaces. 

The initial growth of vapor-deposited copper on differ- 
ent ZnO surfaces at ultra-high vacuum (UHV) conditions 
has been studied experimentally by various techniques 
such as x-ray photoelectron spectroscopy (XPS) 4 i 5 ' 6 ' 7 
low-energy ion-scattering spectroscopy (LEIS) 4 i 5 ' 6 ' 7 
surface x-ray diffraction (SXRD), 8 and scanning tunnel- 
ing microscopy (STM)£ For the two polar surface ter- 
minations of ZnO a growth mode between a Volmer- 
Weber and Stranski-Krastonov behavior was found. Ini- 
tially, the deposited Cu atoms form two-dimensional is- 
lands (i.e. clusters with a thickness of one atom) un- 
til a critical coverage is reached, at which further Cu 
atoms begin to add on top of the islands to form three- 
dimensional clusters, whereas the metal-free regions be- 
tween the clusters fill only very slowlyA^ This behav- 
ior was attributed to purely kinetic limitations at low 
temperatures since a thickening of the islands was ob- 
served upon annealing the surfaces. The critical cover- 
ages at which three-dimensional growth of the clusters 
begins, depends slightly on temperature and was esti- 



mated to be about 1/2 monolayer for the O-terminated 
(0001)-O surface^ and 1/3 of a monolayer for the Zn- 
terminated (OOOl)-Zn polar face of ZnO. 6 In contrast, 
for clean, adsorbate-free nonpolar (1010) surfaces ex- 
clusively three-dimensional islands were observed at all 
coverages £ 

The main focus of theoretical investigations of the 
bonding of metal films and clusters on oxide surfaces 
has been on the MgO (001) surface- 10 ! 11 ! 12 ! 13 ! 14 ! 15 and 
to a lesser extent on the surfaces of more complex ox- 
ides such as Ti0 2f i& AlaOgjiiia and MgAl 2 4 -^ Almost 
all studies were restricted to nonpolar surface termina- 
tions. Polar surfaces, on the other hand, like the two 
basal planes of ZnO, may show a completely different 
behavior compared to nonpolar surfaces. The reason is 
that polar surfaces are characterized by an excess of one 
atomic species 2 ^ which leads to an interesting interplay 
between their electronic and atomic structure: In order 
to be charge-neutral the polar surfaces either have to 
adopt an electronic structure with partially occupied sur- 
face bands (sometimes referred to as "metallization" of 
polar surfaces), or they have to reconstruct (preferen- 
tially by the loss of suitable surface atoms) in order to 
avoid partial band filling, or charged species have to be 
adsorbedi 21 * 2 ^* 2 ^ Depending on which mechanism is real- 
ized in a particular system, very different properties of 
the interaction with metal films and coadsorbates may 
occur. 

Studies of the interaction of copper with ZnO surfaces 
are very scarce. In an early semiempirical quantum- 
chemical studyi the electronic structure and atomic 
charges were calculated for various ZnO clusters with ad- 
sorbed and substituted Cu atom, but no information on 
binding energies and relaxed atomic structures was ob- 
tained. In addition, clusters without embedding were 
used which can not catch the peculiarities of the polar 
surfaces described above. This was taken into account 
in a recent quantum-chemical ab-initio study of the ad- 
sorption of single Cu atoms on different cluster models 24 
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to which we will compare our periodic density-functional 
results. 

In the present study we have calculated the adhesion 
energy and geometry of single Cu atoms and copper 
monolayers at various coverages on the two polar sur- 
face terminations of ZnO. The surfaces are described by 
periodically repeated slabs so that the different stabi- 
lization mechanisms for the polar surfaces, in particular 
the existence of partially occupied surface bands, can be 
handled without any restrictions. Additionally, atomic 
relaxations are fully taken into account. The main focus 
of this study will be on how the charge neutralization 
process of the polar surfaces and the filling of the sur- 
face bands influences the adsorption properties of copper 
spanning the range from single atoms to full monolayers. 
In addition to the very recent quantum-chemical finite 
cluster study, Ref. 0, these are the first ab-initio results 
for the adhesion energies and geometries of copper on 
ZnO surfaces. 



II. COMPUTATIONAL METHOD 

Self-consistent total-energy calculations based on 
density-functional theory (DFT)2£ were carried out to 
determine the adsorption energies and geometries of Cu 
atoms and monolayers on the polar ZnO surfaces. The 
exchange-correlation energy and potential were treated 
within the generalized-gradient approximation (GGA) 
using the functional of Perdew, Burke and Ernzerhof 
(PBE). 26 Normconserving pseudopotentials2£ were em- 
ployed together with a mixed-basis set consisting of plane 
waves and non-overlapping localized orbitals for the O- 
2p, the Zn-3d and the Cu-3<i electrons^ A plane-wave 
cut-off energy of 20 Ry was sufficient to get well con- 
verged results. 

The polar surfaces are represented by periodically re- 
peated slabs. No mirror symmetry parallel to the polar 
surfaces is present which implies that all slabs are in- 
evitably Zn-terminated on one side and O-terminated 
on the other. For the study of Cu monolayers (lxl) sur- 
face unit cells and very thick slabs consisting of 8 Zn-0 
double-layers were used to reduce the residual internal 
electric fields. All atomic configurations were fully re- 
laxed by minimizing the atomic forces. The calculations 
for single Cu atoms were performed on (2x2) surface 
unit cells (corresponding to a Cu surface coverage of 1/4 
monolayer) , and the slab thickness was reduced to 4 ZnO 
double-layers. In a (2x2) surface unit cell arrangement 
the Cu atoms are more than 6.5 A apart. The interac- 
tion energy between the Cu atoms is calculated to be less 
than 0.02 eV so that this setup is sufficient to describe 
very well isolated Cu atoms. 

Three different high-symmetry adsorption sites are 
present at the polar ZnO surfaces (see Fig.^l: an 'on-top' 
position which would be the regular lattice site for the 
next atomic layer, a 'hep-hollow site' above atoms in the 
second surface layer and a 'fcc-hollow site' with no atoms 



(a) top view (b) side view 




FIG. 1: Schematic diagram of the high-symmetry adsorption 
sites on the polar ZnO surfaces, shown here for the oxygen 
surface termination only. 

beneath. All three adsorption sites are considered in the 
present study. Copper atoms are only adsorbed on one 
of the two slab terminations. The electrostatic potential 
is calculated as for a genuinely isolated slab^ 1 '- 1 so no 
artificial interactions between the periodically repeated 
images of the slabs is present. Monkhorst-Pack k-point 
meshes^ with a density of at least (6x6x4) points in the 
Brillouin-zone of the primitive ZnO unit cell are chosen. 
For more details on convergence parameters, the con- 
struction of appropriate supercells as well as the calcu- 
lated bulk and clean surface structures of ZnO we refer 
to Ref. l29l where the same computational settings as in 
the present study were used. 



III. RESULTS AND DISCUSSION 

A. Validation of the Method: Finite Cluster Model 

Before we start our investigation of copper on the po- 
lar ZnO surfaces, we address the question how well our 
DFT method and in particular the PBE functional are 
able to describe the interaction of Cu and ZnO surfaces. 
For Cu on MgO(OOl) significant differences between var- 
ious wave-function based quantum-chemical methods 
and DFT calculations using different GGA functionals 
were observed^ As a general trend it was found that 
DFT slightly overestimates the adhesion energy of Cu 
on MgO clusters compared to best coupled-cluster cal- 
culations, with noticeable variations between the GGA 
functionals. 

To test the reliability of our DFT results we have calcu- 
lated the binding energy and equilibrium distance of Cu 
on a small, isolated Zn404 cluster which we can compare 
with recently published results of an accurate coupled- 
cluster-type calculation (multi-configuration coupled 
electron-pair approximation MC-CEPA)tM A wurtzite- 
like cluster geometry was chosen in order to mimic the 
stacking sequence of polar slabs (see Fig. The Cu 
atom is adsorbed to the top O and bottom Zn atom of 
the cluster which are three-fold coordinated like at a ZnO 
surface. 

The TjU^Oa cluster itself is a closed-shell system, but 
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(a) 





FIG. 2: Adsorption geometries of (a) single Cu atoms and 
(b) Cu dimers on-top of a finite Zn-iCU cluster. 



with the additional Cu atom it becomes open-shell. 
Therefore, in the first calculation for a single Cu atom 
on-top of the three-fold O atom of the cluster we used 
the spin-polarized version of the PBE functional to ob- 
tain the total energy of the adsorbate system as well as for 
the energy of the single Cu reference atom. The binding 
energy and the equilibrium distance are given in Table [I] 
No significant differences in the adsorption properties are 
found if the spin-polarization is completely neglected and 
both, the calculations for the adsorbate system and the 
separate Zn404 cluster and the isolated Cu atom are done 
with the nonpolarized PBE functional. Also using a Cu 
dimer instead of a single Cu atom gave a very similar 
result. Since the inclusion of spin polarization doubles 
the computational cost, we only used a Cu dimer for the 
study of the interaction of Cu with the three-fold coor- 
dinated Zn atom of the finite cluster. 

Compared to MC-CEPA results we find a satisfying 
agreement for the binding energies and equilibrium dis- 
tances. Like in the case of Cu on MgO, the DFT cal- 
culation (here with the PBE functional) slightly over- 
estimates the bond strength. However, the adsorption 
energy is 'equally' to large by 0.15 to 0.25 eV, so that 
the relative stability of the O and the Zn adsorption site 
is quite well described. Overall we find that Cu pref- 
erentially binds to oxygen atoms with an adhesion en- 
ergy of almost 1 eV. The same behavior was also found 
for MgQi 5 -, but with a lower binding energy of 0.4 eV 
(coupled-cluster) and 0.6 eV (DFT-GGA). 



Adsorption 
Geometry 



PBE MC-CEPA 
do [A] £ ad [eV] d [A] £ ad [eV] 



O-Cu 


sp 


1.91 


0.98 


1.95 0.73 




up 


1.91 


0.96 




0-Cu 2 


up 


1.92 


0.92 




Zn-Cu 


sp 






2.80 0.12 


Zn-Cu 2 


np 


2.48 


0.26 





TABLE I: Equilibrium distances do and binding energies 
Sad of Cu atoms and Cu2 dimers on the Zn404 cluster for 
different adsorption sites as indicated; 'sp' refers to spin- 
polarized/open-shell and 'np' to nonpolarized DFT calcula- 
tions. The MC-CEPA values are taken from Ref. 



cluster (see Table and, more importantly, for the non- 
polar ZnO(10l0) surface^. In addition, it is also larger 
than what is typically known from other oxide surfaces. 15 
The reason for this behavior is that at low copper cov- 
erages, the 4s electrons of the Cu atoms fill the partially 
occupied surface bands and thereby stabilize the polar 
surfaces. The same mechanism was found for the ad- 
sorption of atomic hydrogen on the polar O-terminated 
surface (see Ref. H3): at low H coverages the binding en- 
ergy of an H atom relative to a reference of isolated H2 
molecules is more than 2 eV. However, as soon as a H cov- 
erage of 1/2 monolayer is reached, the surface bands are 
filled and the adsorption energy for additional H atoms 
decreases rapidly. As a consequence, surface structures 
with more than 1/2 monolayer of H are not stable in 
thermodynamic equilibrium with a H2 gas phase at any 
temperature and partial pressured Therefore, we can 
expect that the strong binding of copper on the polar 
ZnO surfaces is only possible as long as the surface bands 
are not fully occupied. If the bands are filled, either by 
higher copper coverages or by the presence of coadsor- 
bates like hydrogen, the Cu adsorption properties may 
change drastically. However, for copper there is one sig- 
nificant difference compared to hydrogen: at higher Cu 
coverages strong Cu-Cu bonds will form so that Cu will 
also adsorb at coverages exceeding 1/2 monolayer. 



B. Cu Atoms on Ideal Polar Surfaces 

For a long time it was believed that the polar ZnO 
surfaces exist in an unreconstructed, truncated-bulk-like 
staters If no reconstruction by loss of surface atoms 
and no adsorption of charged species occurs, the polar 
surfaces are inevitably stabilized by partially filled sur- 
face states with a 3/4 and a 1/4 filled band for the O and 
the Zn termination, respectivelyiSiiSSiSiS 

Assuming such a situation in our DFT calculations we 
find a very strong binding of Cu atoms on both polar 
ZnO surfaces (see Table [nj . The Cu binding energies 
are much larger than what is found for the finite Zn^Oi 



C. Cu Atoms on Adsorbate Pre— Covered Polar 
Surfaces 

Several very recent studies have created considerable 
doubt that the polar ZnO surfaces really exist as ideal, 
truncated-bulk-like surfaces. With scanning tunneling 
microscopy (STM) nanoscaled island and pit structures 
were found at the Zn-terminated face^ 5 . which lead to 
an overall loss of roughly 1/4 of the surface Zn atoms. 
Such a stabilization mechanism was confirmed in DFT 
calculations^^ which furthermore predict high concen- 
trations of adsorbed OH-groups in hydrogen-rich envi- 
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(OOOl)-O (OOOl)-Zn 



Adorption Site d [A] £ ad [eV] d [A] E ad [eV] 



on-top 


1.84 


2.57 


2.45 


1.81 


hep-hollow 


1.93 


1.85 


2.08 


1.87 


fcc-hollow 


1.25 


2.87 


1.97 


2.06 



TABLE II: Equilibrium distances do an d binding energies 
Ea_d of single Cu atoms for different adsorption sites on the 
ideal O- and the Zn-terminated polar ZnO surfaces with par- 
tially filled surface bands. Slabs with a (2x2) surface unit cell 
and a thickness of 4 Zn-O double-layers were used for the cal- 
culation. 



ronments. For the O-terminated surface it was shown by 
He atom scattering, 37 in a study of the CO adsorption 
energetics and also by DFT calculations^ 3 - that 
they are H covered for a wide range of temperatures and 
H partial pressures. Both mechanisms, loss of Zn atoms 
and adsorption of hydrogen, lead to a filling of the sur- 
face bands. The filling need not to be perfect, however, 
in order to study the consequences of such a stabilization 
of the polar surfaces on the Cu (co)adsorption proper- 
ties and also in order to support our prediction from the 
last subsection, we have considered two idealized surface 
structures of the polar ZnO surfaces for which the surface 
bands are fully occupied. 

For the O- and the Zn-terminated surface we assume 
that 1/2 monolayer of H atoms and OH groups, respec- 
tively, is present. We place the H-atoms/OH-groups 
on every second atomic row at the surface. At the 
adsorbate-free rows in between these adatoms the Cu 
atoms are placed and the Cu adsorption energies are com- 
puted as in Sec. lHTBl The results are given in Table ITTT1 
We find indeed that the Cu binding energy is signifi- 
cantly reduced compared to the result of the ideal sur- 
face termination with partially occupied ("metallized") 
surface bands. The adsorption energy of roughly 1 eV 
of Cu atoms on the H covered O-terminated surface 
is only slightly larger than what was found for the fi- 
nite closed-shell Zn404 cluster and for the nonpolar ZnO 
surface (having fully occupied bands, i.e. nonmetallic 
character) i24 In addition, a similarly small value was also 
found in a quantum-chemical embedded cluster study of 
the polar O-terminated surface^* in which the embed- 
ding scheme was chosen in such a way that it leads to fully 
occupied surface bands, similar to our hydrogen covered 
surface. 



D. Cu Monolayers on Polar Surfaces 

We turn now to the study of Cu monolayers on the po- 
lar ZnO surfaces. ZnO and bulk Cu have very different 
lattice constants. The equilibrium distance of the surface 
atoms at the polar surfaces is computed to be 3.28 A (ex- 



Adorption Site d [A] E ad [eV] 

(000l)-O surface covered with 1/2 monolayer of H: 
on-top 1.96 1.12 

(OOOl)-Zn surface covered with 1/2 monolayer of OH: 
fcc-hollow 2.50 0.40 

TABLE III: Equilibrium distances do and binding energies 
B a d of single Cu atoms on polar ZnO surfaces with adsorbate 
structures to fill the surface bands (which results in "insulat- 
ing surfaces" for Cu coadsorption) . Same slabs as in Table |n| 
were used. 

perimentally: 3.25 A), and for the Cu bulk lattice con- 
stant we found 3.65 A which corresponds to a nearest- 
neighbor separation between the Cu atoms of 2.58 A (ex- 
perimentally: 3.615 A and 2.56 A, respectively). Epitax- 
ial growth with a (111) orientation of copper would there- 
fore lead to a lattice mismatch of 21 %. Thus, we can not 
expext that expitaxial growth occurs. Experimentally it 
was found that Cu on the polar surfaces is basically un- 
strained, even at very low coverages. 8 Exclusively (111) 
Cu planes were observed with a fixed orientation given 
by the ZnO substrate. 

The realistic treatment of such unstrained copper lay- 
ers on the ZnO surfaces would be computationally very 
demanding since we would have to use prohibitively large 
surface unit cells for our slab. However, in order to get 
a first idea on how the Cu adsorption strength and the 
preferred adsorption sites change if a high concentration 
of copper is present on the polar surfaces we have nev- 
ertheless studied ideal expitaxial Cu (111) monolayers 
on the ZnO surfaces. For such hypothetical films we 
have calculated the equilibrium distance and the work 
of separation which is defined as the energy difference 
between the adsorbed monolayer and separate ZnO slabs 
and metal films. The results are given in Table For 
the O-terminated surface we find that Cu preferentially 
binds on-top of the oxygen atoms. The equilibrium dis- 
tance and the binding energy of 1.23 eV per Cu atom are 
very similar to the results obtained for single atoms on 
the finite Z^Oj cluster and the surface with filled sur- 
face states. Also for the Zn-terminated surface we find 
a good binding of the copper monolayers with a work of 
separation of almost 0.9 eV per Cu atom. This is in agree- 
ment with experiment where a good adhesion of Cu on 
both polar surface terminations was foundi^ In contrast 
to the O-terminated surface, the hollow positions for the 
Cu atoms are more stable than the on-top site. The two 
hollow positions are very similar in energy and may both 
be populated for non-epitaxial films. This may explain 
the experimental observation that the orientation of de- 
posited Cu relative to the substrate is less strict for the 
Zn-terminated face than for the O-terminationiSi^ Im- 
portantly, we also find that the structural relaxation of 
the ZnO substrate can not be neglected. The relaxation 
energy contributes 20 % to 30 % to the binding energy of 
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rigid slab full relaxation 

Adorption Site d [A] E aep [eV] d [A] E scp [eV] 



Adorption Site 



sequence do [A] d [A] E Bcp [eV] 



(0001)-O surface: 



(0001)-O surface: 










on top 


on-top 


1.98 


1.02 


1.94 


1.23 




hep-hollow 


2.21 


0.38 


1.71 


0.66 


(0001)-Zn surface 


fcc-hollow 


1.75 


0.58 


1.36 


0.90 


hep-hollow 


(OOOl)-Zn surface: 












on-top 


2.56 


0.54 


2.50 


0.62 


fcc-hollow 


hep-hollow 


2.13 


0.63 


1.88 


0.84 




fcc-hollow 


2.12 


0.66 


1.89 


0.87 





IlCp 




1 *n 


i 


fee 


1.96 


1.53 


1.50 


hep 


1.86 


1.48 


0.94 


fee 


1.87 


1.48 


0.93 


hep 


1.90 


1.47 


0.91 


fee 


1.93 


1.48 


0.90 



TABLE IV: Equilibrium distances do and work of separation 
E aep (per surface atom) of a Cu monolayer for different ad- 
sorption sites on the polar ZnO surfaces, with and without 
atomic relaxation of the slab. 



these copper films. 



E. Cu Films on Polar Surfaces 

Finally, we briefly study the effect higher copper cover- 
ages than full monolayer on the bonding at the Cu/ZnO 
interface. This will give some insight into the nature of 
the chemical bond at the interface. The higher coverages 
were simulated by adding a second and a third monolayer 
of Cu on-top of the polar surface. For the two and three 
layer films different stacking sequences are possible. As 
shown in Table and Table IVII the energy differences 
between different stacking sequences turn out to be very 
small so that the specific stacking does not play a role in 
the following arguments. 

Overall we find that the same adsorption sites as 
for the monolayer coverage are the most stable ones 
for the thicker Cu films: the on-top position for the 
O-terminated surface and the hollow-sites for the Zn- 
terminated face. For the O-terminated surface the work 
of separation increases with the thickness of the copper 
film, whereas for the Zn-termination a small decrease in 
the binding energy is observed. This behavior as well 
as the preferred adsorption sites can be explained as fol- 
low: For the O-terminated surface the first copper layer 
adopts the role of the missing next Zn layer. Therefore 
the on-top position, which is the next regular lattice site 
of the wurtzite structure, is preferred. Charge is donated 
to the surface O ions to fill the partially occupied surface 
bands, which can more easily occur for thicker copper 
films, and a partially ionic/covalent bond like in bulk 
ZnO is formed. 

On the other hand, for the Zn-terminated surface, the 
topmost Zn layer behaves more like the first layer of 
the adsorbed metal film. Therefore the hollow positions, 
which correspond to an extension of the metal film, are 
the low-energy adsorption sites. The 4s electrons of the 
partially filled Zn surface band contribute to the metal- 



TABLE V: Equilibrium distances do and work of separation 
-E S ep (per surface atom) of Cu double-layers for different ad- 
sorption sites on the polar ZnO surfaces and different stacking 
sequences. 



Adorption Site sequence do [A] £ ad [eV] 



(0001)-O surface: 
on-top 

(0001)-Zn surface: 
hep-hollow 



hep 
fee 

hep 
fee 



1.96 
1.96 

1.87 
1.88 



1.33 
1.32 

0.76 
0.76 



TABLE VI: Equilibrium distances do and work of separation 
E aep (per surface atom) of Cu triple-layers on the polar ZnO 
surfaces for different stacking sequences. 



lie state of the metal film so that the Zn-Cu bond is 
more metallic in character. However, the "integration" 
of the topmost Zn layer becomes weaker for thicker cop- 
per films. 



IV. SUMMARY AND CONCLUSIONS 

Using DFT calculations and periodically repeated 
slabs to represent the two polar ZnO surfaces, i.e. the 
(OOOT)-O and (0001)-Zn faces, we find a pronounced dif- 
ference in the adsorption energy of Cu atoms depending 
on the coverage. Crucial to the proposed understanding 
of this phenomenon is the fact that the ideal, unrecon- 
structed O- and Zn-terminated surfaces feature partially 
filled surface states with a 3/4 and a 1/4 filled band, 
respectively (often referred to as "metallization" of the 
polar surfaces). Up to 1/2 monolayer coverage the ad- 
hesion of copper on such polar surfaces is found to be 
very strong since the 4s electrons of the Cu atoms can 
fill these partially occupied bands and thereby stabilize 
the polar surface itself. However, as soon as these bands 
are filled due to copper (or due to coadsorbates!) already 
present on the surface, the binding of additional Cu to 
the surface becomes weaker and less favorable in com- 
parison to the cohesion energy of Cu in small Cu clusters 
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and thin Cu films. 

Considering the experimentally known coverage- 
dependent change of the growth mode of copper on the 
polar ZnO surfaces (i.e. formation of two-dimensional 
islands until critical coverages of about 50% and 30% 
for the O- and the Zn-terminated surface, respectively, 
is reached at which three-dimensional growth of Cu clus- 
ters sets in) it may very well be that this sudden change 
is not caused by kinetic limitations at low temperatures 
alone. Instead, it might also be driven thermodynamically 
by the significant change in the adsorption energetics of 
copper on such polar surfaces due to successive band fill- 
ing, which in turn is directly coupled to increasing the 
surface coverage. The lower critical coverage found ex- 
perimentally for the Zn-terminated surface may then be 
explained quite naturally by the lower adsorption energy 
of copper found for this particular surface. Additionally, 
the presence of coadsorbates during or prior to copper de- 
position may strongly influence both the critical coverage 



and the morphology of the grown Cu particles. 

Thus, it becomes evident that the polar ZnO surfaces 
allow for a wealth of quite different adsorption scenarios 
depending on parameters such as coverage, defects, coad- 
sorbates etc. Since these parameters typically change in 
complex ways during catalytic cycles one is tempted to 
speculate that the polar character might be at the very 
heart of the catalytic activity of ZnO. 
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